Convergent expansions for properties of the Heisenberg model for CaV 4 



M. P. Gelfand, 1 ! Zheng Weihong, 2 ! Rajiv R. P. Singh, 3 ! J. Oitmaa, 2 and C. J. Hamer 2 

1 Department of Physics, Colorado State University, Fort Collins, Colorado 80523 
2 School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia 
3 Department of Physics, University of California, Davis, California 95616. 



> 
in 

(N 

o 

m 
O 
vo 

ON 
-I— > 



c 

o 
o 



X 



We have carried out a wide range of calculations for 
the S = 1/2 Heisenberg model with nearest- and second- 
neighbor interactions on a two-dimensional lattice which de- 
scribes the geometry of the vanadium ions in the spin-gap sys- 
tem CaVztOg. The methods used were convergent high-order 
perturbation expansions ( "Ising" and "Plaquette" expansions 
at T = 0, as well as high-temperature expansions) for quanti- 
ties such as the uniform susceptibility, sublattice magnetiza- 
tion, and triplet elementary excitation spectrum. Comparison 
with the data for CaVziOg indicates that its magnetic proper- 
ties are well described by nearest-neighbor exchange of about 
200K in conjunction with second-neighbor exchange of about 
100K. 

PACS numbers: 75.10.-b., 75.10.Jm, 75.40. Gb 



Since the discovery of high temperature superconduc- 
tivity in the cuprates there has been much interest in 
two-dimensional antiferromagnetism. From a theoret- 
ical point of view, there has been an extensive search 
for the so-called spin liquid ground state in models with 
realistic interactions. A possible link between such a 
ground state and high temperature superconductivity 
has often been suggested though not demonstrated. For 
square lattice Heisenberg models it is believed that suf- 
ficient further-neighbor interactions will destabilize the 
Neel state and lead to a spin-gap phase; however, such 
a phase typically has spin-Pcierls order and hence is not 
a true spin liquid A spin-gap phase without other 
spontaneously broken symmetries has been found exper- 
imentally in quasi-one-dimensional spin ladders Q and 
theoretically in strongly coupled two-plane models || . 

In light of these developments, it is interesting that 
the quasi-two-dimensional Heisenberg system CaV^Og 
has been found to exhibit a spin gap Q]. Within each 
layer, the spins form a one-fifth depleted square lattice, 
which we denote the CAVO lattice (see Fig. [l]) . Recently 
this system has attracted much attention from theorists 
|^| |l0|. It is now believed that the spin degrees of free- 
dom in CaV40g are in the "plaquette phase," that is, 
there would be no phase transition if the intact squares 
(plaquettes) of the CAVO lattice were adiabatically de- 
coupled. It has been argued Jj| that the superexchange 
between these spins is mediated by oxygen atoms, which 
are situated out of the plane at the center of the squares 
in the CAVO lattice, and that this may lead to significant 
second neighbor interactions: as a first approximation, 
one expects J\ = J[ = 2J 2 = 2J 2 . Furthermore, it has 



been argued that due to its special structure this lattice 
is susceptible to symmetry-preserving lattice distortions 
fiof . Hence it is appropriate to study a Heisenberg model 
with interactions 

H = J\ ^ Si • Sj + J[ ^ S,; • Sfe 

{i,rn) 

where the first sum runs over nearest-neighbor bonds 
within plaquettes, the second over nearest-neighbor 
bonds between plaquettes, the third over second- 
neighbor bonds within plaquettes, and the last over 
second-neighbor bonds between plaquettes. 

We have carried out a variety of high-order conver- 
gent perturbation expansions for this Hamiltonian that 
are comparable in accuracy to the best available quan- 
tum Monte Carlo calculations (which have been limited 
to J2 = J' 2 = 0) and that provide the only reliable re- 
sults for the parameter values most relevant to CaX^Og. 
An outline of the calculations and results is as follows: 
(1) Ising expansion estimates of the sublattice magneti- 
zation and phase boundary for models with J 2 = J 2 = 
are consistent with the quantum Monte Carlo calcula- 
tions of Troyer et al. ||; (2) Plaquette expansions for the 
triplet elementary excitation spectrum are highly con- 
vergent and informative in the plaquette phase, and (3) 
High-temperature expansions for the susceptibility, com- 
bined with plaquette expansion results, can be compared 
with the existing experimental data Q] to fix the model's 
parameters for CaV40g. In this Letter we can only de- 
scribe the highlights of the results; details are reserved 
for a separate publication. 

For the Ising expansions we have considered the Hamil- 
tonian 
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where e, = ±1 on the two sublattices. (The last term is 
a local field term, which can be included to improve con- 
vergence and proves useful in that regard when J[k. J\.) 
We have developed expansions up to order a 14 , for sev- 
eral values of J[/Ji, for the ground state energy, the 
sublattice magnetization M, and the longitudinal and 
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transverse susceptibilities. The expansion techniques 
have been reviewed earlier j[l],[l2). When the Heisen- 
berg model has a spin gap, one expects these expansions 
to break down at some a < 1 and the system to exhibit 
d = 3 Ising criticality. Such considerations can be used 
to estimate the domain of long-range Neel order. Within 
the Neel-ordered phase our estimates of M are presented 
in Fig. ^ along with those of Troyer et al. ||. 

For the plaquette expansions we return to the orig- 
inal Hamiltonian of interest, Eq. (|l|), but with the in- 
terplaquette interactions parametrized as J[ = \\J\, 
J 2 = 7J1, and J 2 = X 2 J 2 . For convenience we have 
taken Ai = X 2 = A in our calculations; that restriction 
could be relaxed if desired. Let us restrict our attention 
to 7 < 1, since there is both a ground state and low- 
est excited state level crossing for an isolated plaquette 
when J\ = J 2 , and J\ > J 2 is certainly more relevant 
to CaV40g. When A = the ground state is a prod- 
uct of singlets on the non- interacting plaquettes, while 
the lowest-lying triplet excitations are gapped and com- 
pletely local. For A ^ 0, the excitations are mobile and 
develop a well defined dispersion. Using a recently devel- 
oped method for computing the excitation spectra [ p^[ , 
the complete dispersion for these triplet excitations has 
been constructed to order A 5 for various values of 7. Fol- 
lowing the convention of Katoh and Imada [0] for the 
wave vectors, to second order in A the dispersion has the 
form 

A(^, % )/J 1 = l + A 2 c f 0) (7) 
+ \\(cosq x +cosq y )[c[ 1 '°\j) + \c£' 0) ("/)} 

, \2 (1,1)/ \ 

+A cos q x cos q y c\ '(7) 
+ iA 2 (cos2 fe + cos2 9 j,)4 2 ' 0) ( 7 ) (3) 

The coefficients c r n for 7 = 0, 0.4 and 0.5 are presented 
in Table Q. (Because 7 is a parameter in both the un- 
perturbed and perturbing Hamiltonians, a simple double 
expansion in which c r n can be expressed as a polynomial 
of order n in 7 does not exist.) Second-order excitation- 
spectrum expansions have been presented earlier: for 
J 2 = by Katoh and Imada Q and at q x = q y = tt 
for J 2 — J 2 (that is, Ai and 7 arbitrary and A2 = 1) 
by Ueda et al. ||. However, in both cases the authors 
failed to consider the 0(A 2 ) term associated with the 
(cosq x + cosq y ) part of the dispersion; in effect they set 

Cj 1,0 "* to zero. For the case A = 1 and 7 = one finds 
from summing the correct series terms to 0(A 2 ) a gap 
A(7r,7r) = O.OlOJi, in contrast to the previous pertur- 
bative estimates of 0.205Ji. As evidenced by the Ising 
expansions and the Monte Carlo results of Troyer et ai, 
this point lies in the Neel-ordered phase and should have 
a zero gap. This shows that already in second order the 
plaquette expansion for the gap locates the phase bound- 
ary fairly well, at least for 7 » 0. 

Another important point not discussed in earlier 



perturbation-theoretic work is that the minimum of A 
is not necessarily found at q = (tt,tt). For small A the 
minimum shifts from (tt, it) when 7 < 1/2 to (0,0) when 
7 > 1/2. Precisely at 7 = 1/2 the minimum for small 
A is achieved at (q, q) with q/n = 0.7044 .... To the ex- 
tent one can trust estimates of A based on simply sum- 
ming series coefficients, such an incommensurate mini- 
mum persists at finite A and for 7 w 1/2. (A similar 
incommensuration for sufficiently large 7 was found in 
the mean-field calculations of Starykh et al. JlO|.) 

The present plaquette expansions are also the first 
which are sufficiently long to reflect the absence of re- 
flection symmetry in the Hamiltonian: in third order and 
above A(q x , q y ) =^ A(g a , q x ). It is not a dominant feature 
of the excitation spectrum but if 7 is sufficiently large it 
may be observable: at 7 = 0.5 and A = 1 the asymmetry 
parameter 

2[A(g, tt - q) - A(tt - q, q)]/[A(q, tt - q) + A(tt - q, q)] 

(4) 

may be 10% for q/ir w 0.15. 

The fact that the excitation spectrum minimum is not 
fixed at q = (tt, tt) could in principle complicate the mat- 
ter of estimating the gap in the plaquette phase. How- 
ever, for 7 < 1/2 and A < 1 the value of A(tt, tt) exceeds 
the true gap by at most 0.05 Ji and usually by far less. So 
it is not too bad to identify estimates of A(tt, tt) with esti- 
mates of the gap in the parameter region of most physical 
interest. Another consideration is that the uncertainty in 
the estimates of A(tt, tt), which can be judged by spread 
in the direct sums and Pade approximants, become no- 
ticeable precisely in the part of 7-A space where the gap 
minimum apparently moves away from q = (tt, tt) . See 
Fig. H for some examples of A(tt, tt) estimates. Fortu- 
nately, for the parameters expected to be relevant to 
CaV40g the estimates of the gap appear to be fairly re- 
liable, if not quite as good as when 7 r* 0. For those 
parameters, J[ — J\ — 2J 2 = 2J 2 (A = 1, 7 = 1/2), the 
gap is 0.57±0.03 J\. Taking the experimental gap value of 
107 K for CaV 4 9 leads to the estimate J x = 190±10K. 
The Curie temperature (= 9Ji/8, for the given parame- 
ter ratios) should then be 215 ± 10 K, which is consistent 
with the experimental value of 220 K and shows that the 
present parameter estimates are entirely plausible. 

If magnetic neutron scattering data can be obtained, 
comparison with the plaquette expansion for the gap 
should allow for precise determination of model param- 
eters. The structure of the triplet excitation spectrum 
is strongly dependent on the ratio of nearest- to second- 
neighbor couplings: see Fig. ^ for some illustrative plots. 

We now turn to the high temperature expansions for 
the uniform susceptibility. The uniform susceptibility 
per spin can be defined as /3x(T) — ^J2iJ2j(^i^j) 
where the angular brackets represent thermal averaging 
and (3 = 1/T. For J 2 = J 2 = we have determined 
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the expansion to order /3 14 . The integrated differential 
approximants have been used to extrapolate \ to low T; 
we find that they agree with Monte Carlo calculations 
Q down to T » Ji/10 when the ground state lies in the 
plaquette phase. 

High temperature expansions have also been developed 
for the models with J\ = J[, J 2 = J 2 i= 0, to order 
(3 W . They should be more relevant for comparison with 
the CaV40g data, however, they have proven more dif- 
ficult to extrapolate to lower temperatures. In order to 
improve the convergence, we have developed differential 
approximants for xe A / T , with A taken from our studies 
at T = 0. This biases the leading exponential behavior of 
the susceptibility at low temperatures to be of the form 
g— a/t^ D eS pit e this biasing there remain significant un- 
certainties in the location of the peak temperature, where 
X is maximized; the best we can say is that it lies in the 
range 0.4-0.6 J\ for J 2 = Ji/2. The susceptibility of 
CaV40g peaks at about 110 K, which is consistent with 
a Ji value of approximately 200 K. 

A comparison of experimental data with the calculated 
susceptibility, using J\ — 190 K as estimated above, is 
shown in Fig. ^. In order to convert the susceptibility 
for the Heisenberg model into units of emu/g, we mul- 
tiply the theoretical value (calculated with J\ = 1) by 
(.9 2 'Mb I '</i&B)(4iV a/ 'M) , where /z_b is the Bohr magneton, 
Icb is the Boltzmann constant, Na is Avogadro's number, 
M is gram molecular weight, and the factor of 4 comes 
from the four vanadium spins per CaV^g formula unit. 
The <?-factor can be used as a fitting parameter, and the 
comparison shown uses the value 1.77. We note that 
the calculated susceptibility does not fully agree with the 
data at high temperatures: this may be in part due to 
the existence of other phases above 340-ftT as noted in 
the experimental paper W, or due to the lattice distor- 
tions discussed by Starykh et al. (To), which can make 
the exchange constants temperature dependent. A more 
detailed examination of the high temperature suscepti- 
bility for CaV40g in context of the Heisenberg models 
will be given elsewhere. 

In summary, we have carried out a wide range of high- 
order perturbation expansions for S = 1/2 Heisenberg 
antiferromagnets on the CAVO lattice. The existing sus- 
ceptibility data for CaV40g is generally consistent with 
the coupling ratios J\ = J[ = 2J 2 — 2J' 2 suggested by 
Ueda et al. (|) . Measurements of the wave vector resolved 
triplet excitation spectrum should be compared with es- 
timates of the spectrum based on the present plaquette 
expansions for a definitive test of the model. 
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FIG. 1. The CAVO lattice, with sites indicated by cir- 
cles. The couplings J\, J[, J2, and J' 2 are indicated by thick 
solid, thick dashed, thin solid, and thin dashed lines, respec- 
tively. Note that the plaquette centers lie on a square lattice 
with spacing b which is y/Z times the distance between near- 
est-neighbor sites. In characterizing the excitation spectrum 
(see Eq. ([|) below) we take 6 = 1 and rotate the coordinate 
system so that the lines between nearest-neighbor plaquette 
centers define the x and y axes. 



FIG. 2. Sublattice magnetization versus J'i/Ji (with sec- 
ond-neighbor couplings J2 = J2 = 0) as estimated by Ising 
expansions (filled symbols) and the Troyer et al. Q quantum 
Monte Carlo calculations (open symbols). 



FIG. 3. The excitation energy A(n,ir) versus A for several 
values of 7. 



FIG. 4. The triplet excitation spectrum A(q x ,q y ) along 
high-symmetry cuts through the Brillouin zone for coupling 
ratios A = 1 and 7 = 0.1, 0.3 and 0.5 (from greatest to least 
values of A(0, 0)). The ends of the bars indicate the values of 
the direct sums to fourth and fifth order in A. 
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FIG. 5. A comparison of the calculated susceptibility for 
the Heisenberg model with J2 = Ji/2 (the solid lines repre- 
senting various approximants), with the experimental data of 
Taniguchi et al. j| for CaV^Og. 

TABLE I. Coefficients in the expansion of the triplet ele- 
mentary excitation spectrum A(q x ,q y ), to second order in A, 
for selected 7 values: see Eq. (J3J). 



7 = 



7 = 0.4 



7 = 0.5 




0.0358796 
0.6666666 
0.1944444 



0.0185185 



0.1111111 



0.0444444 



0.1950309 
0.1333333 
0.1799242 



0.0570370 



0.1597222 
0. 



0.2279762 
0. 



0.0666667 
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